Light wheel buildup using a backward surface mode 
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When a guided mode is excited in a dielectric slab coupled to a backward surface wave at the interface between 
a dielectric and a left-handed medium, light is confined in the structure : this is a light wheel. Complex plane 
analysis of the dispersion relation and coupled-mode formalism give a deep insight into the physics of this 
phenomenon (lateral confinement and the presence of a dark zone). 
OCIS codes: 260.2110, 160.4670. 



Recently, a lamellar structure consisting of a conven- 
tional dielectric layer coupled to a left-handed material 
(LHM) layer has been proposed to confine light. Because 
of contradirectional power flows in the two layers, an 
exotic localized mode called "light wheel" can be ex- 
cited [1,2]. In this letter, we present our study of a new 
type of light wheel that uses a backward surface mode at 
the interface between an LHM and a right-handed mate- 
rial (RHM) . The dispersion relation of the structure and 
the coupled-mode formalism allows us to describe and 
explain the field distribution in the structure. 

Let us first consider a surface wave propagating along 
an interface between two semi-infinite media: an LHM 
and a RHM. The LHM is caracterized by its dielectric 
permittivity £3 and its magnetic permeability [13, which 
are both negative. The permittivity of the RHM is £1, 
its permeability is fj,\. Because of unusual properties of 
LHM, such an interface supports surface guided modes 
which can be backward (i.e., present opposite phase and 
group velocities) [3,4]. In TM polarization, the disper- 
sion relation for surface waves propagating along this 
interface can be written as [4] 
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where X = J^i, 

and fco is the wavenumber in the vacuum. A backward 
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^0-, a is the propagation constant 



surface wave is supported provided the inequality 
1 < X < 1/Y 



(2) 



is verified. 

The light wheel develops from the contradirectional 
coupling between this backward mode and a forward 
guided mode. This coupling appears when the disper- 
sion curves of the two waveguides cross each other (i.e. 
when they are under phase matching conditions). Even 
when taking the dispersive character of the LHM into 
account, such a coupling is thus not difficult to obtain 
as long as the condition (2) is fulfilled. At a given wave- 
length A, let us consider, for instance, that £3 = —0.8 
and /A3 = —1.5. Because that satisfies the previous con- 
ditions, a backward surface wave between such an LHM 



and the air (£1 = 1, /ii = 1) exists for a = ao = 1.2472fcg 
according to Eq. (1). 

Let us obtain the coupling with a symmetrical dielec- 
tric slab waveguide as depicted Fig. 1. According to its 
dispersion relation, the thickness of the waveguide sup- 
porting a forward guided mode for the propagation con- 
stant ao, in TM polarization, is given by 
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; — £1 fJ-i fc„ and 72 = \J £2 A*2 fco - a\ 



where k% = \/ 
For E\ = 1, fj,i = 1, £2 = 3 and /x 2 = 1, the value of 
thickness /i 2 is 0.2854A. In the following, we retain the 
above characteristics. 

The above dielectric slab and LHM interface arc 
brought close together, forming the structure presented 
in Fig. 1. We have calculated the analytical expression 
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Fig. 1. Dielectric slab waveguide and the LHM interface 
separated by a distance hi and surrounded by a medium 
characterized by £1 and /xi. 

of the dispersion relation of the whole structure, in TM 
polarization: 
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[cxp(2j7 2 fo 2 ) + X 2 ] [cxp(2wi/n) + A 3 ] 
[exp(2j7 2 /i 2 ) - A 2 ][exp(2Ki/ii) - X 3 ] 
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where Xi = KiE\j K\Ei and A^ = (1 — Xi)/(1 + Xi). Figure 
2 shows, in the alpha complex plane, the solutions of the 
dispersion relation (Eq. (4)) for several values of distance 
hi, the frequency being fixed. When distance hi between 
the slab and the interface is large enough, the waveguides 
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are independent, the dispersion relation is thus verified 
for a = olq. But when distance hi decreases, they become 
coupled and two complex-conjugate solutions appear. 
For instance at h x = 0.8A, a = (1.247994±0.012785i)fc . 
These twin solutions have the same real part so that the 
corresponding modes cannot be excited separately by a 
source and they form a light wheel as shown in Fig. 3. 

The nonzero imaginary part of a means that the field 
E(x,z,t) = Eo(z,t) exp(iax) decays along the x axis. 
Because the imaginary parts of the solutions are oppo- 
site, the twin modes decay in opposite directions. The 
size of the light wheel is moreover controled by 9(a), 
a characteristic width of the light wheel being given by 
L = cj^jy- Figure 2(a) shows that the smaller the dis- 
tance hi is, the higher the imaginary part of the propa- 
gation constant is and so the smaller the light wheel will 
be. However, Fig. 2(b) shows that the imaginary part of 
a presents a maximum reached for strong coupling for 
hi = 0.17A. The light wheel can, thus, confine the light 
in a region as small as a characteristic width L = 1.17A. 
Figure 3(a) shows the electromagnetic field created in the 




Fig. 3. Modulus of the field represented in a domain f 00 A 
large, 5 A high. The distance between the dielectric slab 
and the LHM interface is hi = 0.8A. (a) The punctual 
source, placed in the dielectric waveguide, excites two 
contrapropagative light wheels, (b) The light wheel is 
excited by a an incident gaussian beam (angle: 33.9°, 
waist: fOA ) in a prism (e = 5, fi = I). White arrows 
indicate the propagation direction of light. These images 
are obtained using the numerical method described in [9] . 
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Fig. 2. Solutions of the dispersion relation in the a com- 
plex plane, for different values of hi in wavelength, the 
distance between the dielectric waveguide and the LHM 
surface. Inset, zoom in the region around ao- 

structure by a punctual source inside the dielectric layer. 
The interference pattern is due to the two contrarotative 
light wheels, which are excited and interfere. In Fig. 3(b), 
the light wheel is excited by a beam using evanescent 
coupling: a prism is placed above the dielectric slab and 
the structure is illuminated by a gaussian beam coming 
from above. Here, the field distribution is less intuitive. 
In particular, a dark zone appears just below the inci- 
dent beam, in the center of the dielectric waveguide but 
not at the plasmonic interface. The presence of a similar 
dark zone has already been noticed but not explained [1] . 
We have applied the coupled-mode theory (CMT) [5-7] 
to this contradirectional coupling to get an analytical 
model able to account for this phenomenon. 

Let us consider two independent guided modes whose 
complex amplitudes are A and B. Mode a corresponds 



to a right traveling wave in the dielectric slab, whereas 
mode b travels to the left on the LHM interface. If the 
two guides are brought close together, they become cou- 
pled: energy is exchanged between them. Hence in this 
contradirectional coupling case and under phase match- 
ing conditions, the complex amplitudes A and B depend 
on x, obeying relations of the type [6] 

(5) 
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= kA, 



(6) 



where k is the coupling coefficient. Its value is given by 
the imaginary part of the solution a of dispersion relation 
Eq. (4) . For weak coupling as well as for strong coupling, 
we indeed have k = iQ(a). 

Here, the evanescent coupling used to excite mode a 
is assumed to be weak so that the guided modes are 
unmodified by the presence of the prism. An approach 
suggested by [8] is to consider the guided mode excited 
by a set of punctual sources situated inside the waveg- 
uide with an amplitude distribution along the x direction 
given by the incident beam. Each punctual source has an 
amplitude proportional to the amplitude of the incident 
field at the prism interface just above. 

Let us then determine A x (x) and B\(x), the modes 
amplitudes created by a punctual source localized at x = 
whose the amplitude is equal to one, i.e., a source given 
by the expression S±(x) = S(x), where 6 is the Dirac 
distribution. A\ and Bi can be seen as Green's functions 
[8]. The solutions of Eqs. (5) and (6) are 

Ai(x) = H(x) exp(-|«|a;) - H(-x) exp(|/c|a;), (7) 

Bi(x) — exp (^— i^j [H(x) cxp(— \k\x) + H(—x) cxp(|k|x)] 

(8) 
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where H is the Heaviside function. Notice the 7r/2 phase 
shift between the amplitudes A\ and B\ in each domains 
x > and x < 0. 

When the structure is excited by a source amplitude 
distribution S(x), the amplitudes of the modes are ob- 
tained by convolution of the source expression S with the 
Green's functions Ai and B±. Thus, for a gaussian beam 
whose the waist size is 10, i.e., an amplitude distribution 
S(x) = So cxp (— x 2 /w 2 ), we get 
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where erfc is the complementary error function. 

Figure 4 shows the modulus of the field distribution in 
the middle of the dielectric slab of Fig. 3(b). It is com- 
pared to the theoretical mode profile computed 
for k = 0.080330i, given by the dispersion relation. Be- 
cause the coupling with the prism is weak, the coupled 
modes arc relatively undisturbed and Eq. (9) accurately 
describes the field in the slab. Both curves show a min- 
imum for x = which corresponds to the dark zone in 
the middle of the dielectric waveguide. The light wheel 
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Fig. 4. Modulus of the field in the middle of the dielec- 
tric waveguide (blue curve) and |A(a;)| the modulus of 
the theoretical amplitude of the mode a (red curve). Pa- 
rameters of the structure and the beam arc as in Fig. 3. 
So is chosen arbitrarily. 



contradirectional coupling to the plasmonic backward 
mode with a — 7r/2 phase shift. There is no reason why 
the backward mode should undergo any phase change 
in x = 0, and the energy is transferred to the dielectric 
waveguide for x < 0, with another phase shift of — 7r/2. In 
the dielectric slab, the right part and the left part of the 
light wheel are thus in phase opposition. This cannot be 
seen in Fig. 3(a) because the source is ponctual. In Fig. 
3(b), because the evanescent coupling is equivalent to a 
spatially extended source, the parts of the light wheel 
which arc in phase opposition "overlap" in the dielectric 
slab and a dark zone appears. 

When we consider a lossy LHM, which is more likely, 
the light wheel and the dark zone still exist as shown 
in [1] . In this case, the model can be simply extended by 
changing Eq. (6) into dB/dx = kA — kiB, where Ki is 
the extinction coefficient. 

Finally, the complex plane analysis and the CMT as- 
sociated with Ulrich's approach of evanescent coupling 
allow a very accurate description and a deep understand- 
ing of the light wheel phenomenon and its universal fea- 
tures. Beyond the fact that a light wheel can be used 
to confine light, this phenomenon can be used for beam 
reshaping [9], and an analytical model is in this context 
particularly useful. 
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phenomenon can be summarized from Eqs. (9) and (10) 
as follows: the guided mode is excited in the dielectric 
waveguide towards the right. It is then transferred by 
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